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Abstract 

We analyze the anomalies of superconducting state (s and d-wave pairing) in 
a simple model of pseudogap state, induced by fluctuations of short - range 
order (e.g. antiferromagnetic) , based on the model Fermi surface with "hot 
patches". We derive a system of recursion relations for Gorkov's equations 
which take into account all diagrams of perturbation theory for electron inter- 
action with fluctuations of short - range order. Then we find superconducting 
transition temperature and gap behavior for different values of the pseudo- 
gap width and correlation lengths of short - range order fluctuations. In a 
similar approximation we derive the Ginzburg ~ Landau expansion and study 
the main physical characteristics of a superconductor close to the transition 
temperature, both as functions of the pseudogap width and correlation length 
of fluctuations. Results obtained are in qualitative agreement with a number 
of experiments on underdoped HTSC - cuprates. 

PACS numbers: 74.20.Fg, 74.20.De 
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I. INTRODUCTION. 



The pseudogap state, observed mainly in the underdoped region of the phase diagram of 
high - temperature superconducting (HTSC) cuprates, leads to a wide range of anomalies of 
physical properties of these system both in normal and superconducting states |l|]. There are 
two main theoretical scenarios to explain these anomalies. The first is based upon the model 
of Cooper pairs formation already above the temperature of superconducting transition 
with phase coherence appearing only for T < T^. The second assumes that the pseudogap 
state is induced by fluctuations of antiferromagnetic (AFM) short - range order which exist 
in the underdoped region of the phase diagram 10-0. Recently a number of experiments 
appeared which rather convincingly favor the second scenario 

Most theoretical papers in the field are dedicated to the study of the models of the 
pseudogap state in normal phase for T > T^. In Refs. |]Tn|,|n[ we proposed a very simple 
(toy) exactly solvable model of the pseudogap state, which assumes the existence of certain 
"hot" (flat) patches on the Fermi surface. Within this model we derived the Ginzburg - 
Landau expansion for different types of Cooper pairing [|10| and studied the anomalies of 
superconducting state for T < Tc induced by fluctuations of AFM short - range order. 
In these papers we dealt with (over) simplified model of static Gaussian fluctuations of short 
- range order with infinite correlation length, which allowed us to obtain an exact analytical 
solution for the pseudogap state. In real systems the correlation length of AFM fluctuations 
is finite and relatively short 0. The main aim of the current paper is to generalize the 
main results of Refs. [0,|n|] to the case of finite correlation lengths of fluctuations of AFM 
short - range order and to analyze the dependence of the main physical characteristics of 
the superconducting state on this correlation length, as well as on the effective width of the 
pseudogap. 



II. MODEL OF THE PSEUDOGAP STATE. 



We consider the simplified model of the pseudogap [|10|,|rTl based upon the picture of 
well - developed fluctuations of AFM short - range order and is in some respects close to 
the "hot - spot" model of Ref. 0. We assume the Fermi surface of two - dimensional 
electronic system as shown in Fig.l. In fact such a Fermi surface was observed in a number 



of ARPES - experiments on HTSC - cuprates [|T^,|T^. Note that the assumption of the 
flatness of these patches is non crucial for our model, but significantly simplify calculations, 
which, in principle, can be done also in more realistic "hot-spot" approach with apparently 
similar results. The model Fermi surface as in Fig.l has already been applied to HTSC 
- cuprates in Refs. |]T3|-|I^, where the details of microscopic criteria for the existence of 
antiferromagnetic and superconducting phases were analyzed. Here we just assume a kind 
of phenomenological model with static Gaussian fluctuations of short - range order with 
correlation function (structure factor) of the following form 0: 

where ^ - is correlation length of fluctuations, while the scattering vector is taken to be 
Qx = ±2kF, Qy = or Qy = ±2kF, = 0, anticipating that these fluctuations are 
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incommensurate. The factorized form of (||) was first introduced in Ref. and leads 
to great simplification of all calculations. At the same time we stress that in practice it 
practically coincides with the usual isotropic Lorentzian in the most important region of 

|q-QI <r' 0- 

Less justified physically is the assumption of the static nature of fluctuations which can 
be reasonable only at high enough temperatures [Pfl- temperatures, including the 

superconducting phase, spin dynamics can be very important, e.g. for the microscopies of 
Cooper pairing according to the picture of "nearly - antiferromagnetic" Fermi liquid p!7| , p!8| . 
However, we assume that our static approximation can be sufficient for qualitative analysis 
of the influence of superconductivity upon superconductivity, which will be described below 
within the standard BCS - like approximation. 

Let us write down the interaction of electrons with AFM fluctuations in the following 
form: 

Veff = {27^fW^S{ci) (2) 

where W is determining the energy scale (width) of the pseudogap. We assume that only 
electrons from the flat ("hot") of the Fermi surface interact with AFM fluctuations, so 
that W is effectively non zero only for these electrons pO|JTl[] . Note that we completely 
neglect the spin structure of interaction which can be rather easily taken into account 0, 
but makes calculations more cumbersome. In this sense, strictly speaking our discussion is 
more appropriate for the case of electrons interacting with fluctuations short - range order 
of CDW rather than SDW (AFM) type. We also assume that this simplification is relatively 
unimportant for the analysis of qualitative effects of pseudogap on superconductivity. 

The factorized form of the correlator (|lD and of interaction (0) leads to the one - dimen- 
sional nature of scattering on fiuctuations. In the limit of infinite correlation length ^ ^ oo 



this model acquires an exact solution |T0| , p!T|J19|| . For finite ^ we can construct "nearly exact" 



solution [|^, directly generalizing the one - dimensional approach proposed in Ref. |2^. In 
this case we can (approximately) sum the whole diagram series for one - particle Green's 
function of electrons from the fiat parts of the Fermi surface (where the "nesting" condition 
for the electronic spectrum ^p±Q = — is satisfied). 

For the contribution of an arbitrary diagram for the self - energy of an electron of the 
A^-th order in (Q) we use the following Ansatz |17|,|20[|: 



2N-1 



=1 



where k = vfC,~^ {vp - Fermi velocity), kj - is the number of interaction lines surrounding 
the j-th (from the beginning) electronic line in diagram, e„ = 27rT(n+ 1/2) (in the following 
we write expressions for the case of > 0). Thus the contribution of an arbitrary diagram 
is determined, in fact, only by some set of integer numbers kj. An arbitrary diagram with 
intersections of interaction lines is seen to be equal to some diagram without intersections 
and the contribution of all diagrams with intersection can be taken into account with the 
help of some combinatorial factors v{kj), attributed to interaction lines on diagrams with no 
intersections In the model of incommensurate fiuctuations which we only consider: 
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v(k) = I ^ ^"^"^ ^ (4) 
1 I for even k ^ ' 

As a result we obtain the following recursion procedure (continuous fraction representation) 
to calculate the one - particle Green's function G(e„p) of electrons from "hot patches" 



G^-(-»P) = _ + - WMk + l)G.«(.„p)' G(.„p) . G„(.„p) (5) 

Diagrammatically this procedure is shown in Fig.2. 

Ansatz (|^) for the contribution of an arbitrary diagram of the A^-th order, in general 
case, is not exact ||7,21]. However, in two - dimensional case we can explicitly demonstrate 



the topologies of the Fermi surface when becomes exact |]^, for the general case we can 
show 1^ that this representation in some sense overestimates the role of finiteness of corre- 
lation length ^ in the given order of perturbation theory. For one - dimensional case, when 
this problem is particularly serious |17|,|2T[|, it appears that for incommensurate fluctuations 



calculations of the density of states using @ quantitatively almost ideally reproduce |22 
the results of an exact numerical solution given in Refs. p3| , |2^ [|. In the limit of ^ cxo 
Ansatz (§) reduces to an exact solution of Ref. [|l9l, while in the limit of ^ — » for the fixed 
value of W it describes the physically correct limit of free electrons. 

Outside "hot patches" electrons in our model do not interact with fluctuations and the 
Green's function remains a free one: 

G{enp) = GooienP) = -. — ^j— (6) 

This model leads to non Fermi - liquid like behavior of the spectral density on "hot patches" 
and to the smeared pseudogap in the density of states (cf. similar results in "hot spots" 
model P,^). On "cold" parts of the Fermi surface we just get the usual Fermi - hquid (free 
electron) behavior. 



III. GORKOV EQUATIONS FOR A SUPERCONDUCTOR WITH PSEUDOGAR 



In Refs. 10.11 



we have analyzed the anomalies of superconducting state in an exactly 
solvable model of the pseudogap state induced by AFM short - range order fluctuations with 
infinite correlation length (^ oo). In particular, in Ref. ||lT| it was shown that AFM fluc- 
tuations may lead to strong fluctuations of superconducting order parameter (energy gap A) 
breaking the standard assumption of self - averaging gap p5|-p7| which allows independent 



^In the case of one - dimensional problem with commensurate fluctuations Ansatz <^ does not 
describe a certain weak Dyson - type singularity of the density of states near the center of the 
pseudogap |23,24|, while outside the region of this singularity it also gives rather good quantitative 
description of exact results. Note that in two -dimensional case Dyson - singularity, most probably, 
is just absent. 
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averaging (over configurations of tlie random field of static sliort - range order fluctuations) 
of the order parameter A and different combinations of electronic Green's functions, enter- 
ing the basic equations of the theory. Usually, the possibility of such independent averaging 
is supported by the following argument [^,^: the value of A changes on characteristic 
scale of the order of superconducting coherence length ,^0 ~ vf/Aq of BCS - theory, while 
Green's functions typically vary on much shorter lengths of the order of interatomic dis- 
tances. Naturally, the last assumption becomes wrong if a new length scale ^ — > cxd appears 
in electronic system. However, in case of AFM correlation length ^ <C ^0 (i-e- when AFM 
fluctuations correlate on distances much smaller than the size of Cooper pairs) the self - 
averaging property of A remains, being broken only for C, > ■Co- For this reason all our 
analysis below will be done assuming the self - averaging nature of the energy gap, allowing 
us to use the standard approach of the theory of disordered superconductors (mean - field 
approach in terms of Ref. [ll|). Thus, we temporarily put aside the very interesting problem 



of superconductivity in the absence of self - averageness of the order parameter. Note that 
in real HTSC - systems we, in fact, usually have ^ ~ so that these materials are in the 
region of parameters most difficult for the theory. 

Similarly to Refs. we assume that superconducting pairing is due to the following 

simplest attraction potential: 

V{p,p') = V{<P,<j>') = ~Ve{<jy)e{<f>'), (7) 

where (p is the polar angle determining the direction of electronic momentum p in highly 
conducting plane, while for e(0) we take the simple model dependence: 

, , J 1 ( s-wave pairing) , . 

I v^cos(20) ( d-wave pairing) ' 

The attraction constant V is, as usual, non zero in some shell of the width of 2ujc around the 
Fermi level [ujc ^ is characteristic frequency of the quanta responsible for electron attraction). 
In this case superconducting gap takes the form: 

A(p)=A(0) = Ae(0). (9) 

In the following for brevity we shall denote gap simply as A instead of A(0) writing down 
explicit angular dependence only where it is necessary. 

In the superconducting state perturbation theory over the interaction with AFM fluctu- 
ations (IID has to be built on "free" normal and anomalous Green's functions of a supercon- 
ductor: 

In the model with flat patches on the Fermi surface the electronic spectrum on the patches 
orthogonal to is: = f^dPsI — Pf) as electron velocity v is perpendicular to Py (every- 
thing is symmetric for patches orthogonal to Py). Thus, in the case of s-wave pairing, when 
A is independent of the direction of electronic momentum, for the model interaction (|l|), (^ 
the problem becomes one - dimensional. In case of d-wave pairing situation is more difficult 
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as even on the flat patches orthogonal to px the value of A(0) depends on py (and symmet- 
rically on patches orthogonal to Py). For this reason, in the analysis of d-wave pairing it is 
convenient to introduce instead of (0) the fluctuation correlator of the form 

S{cL) = -\ , ^/ L + , ^/ !2^.-2 ^(g-) I (11) 

In this case interaction does not change Py or p^ on flat patches orthogonal to px or py and 
the problem again reduces to one - dimensional. 

Now we can formulate an analogue of approximation (^) also for superconducting state. 
Some details of derivation are given in the Appendix A. The contribution of an arbitrary 
diagram of A^-th order over interaction (0) to normal or anomalous Green's function is 
given by the product of A^ + 1 "free" normal Gofcj or anomalous Fqj.. Green's functions with 
renormalized frequencies and gaps (cf. below). Here kj - is the number of interaction lines 
surrounding the j-th (from the beginning of the diagram) electronic line. As in normal phase, 
the contribution of an arbitrary diagram is defined by the set of integers kj and each diagram 
with intersection of interacting lines is equal to a certain diagram of the same order without 
intersections. Thus again we can use only diagrams without intersections of interaction lines 
accounting diagrams with intersections by the same combinatorial factors v{k) (attributed 
to interaction lines) as in the normal phase. As a result we obtain diagrammatic analogue 



of Gorkov equations |28| shown in Fig. 3. Accordingly we get a system of two recursion 



relations for normal and anomalous Green's functions: 

Gk = Gok + GokGGk — GquFF^ — FQkG*F^ — F^^F^Gk 



k ^ Ok^ ^ k 

Fk = Fok + FokGGk - Ff^^F F^ + Gq,,G*F^ + Gl^^F^Gk (12) 



where 



G = W\{k + l)Gfc+i; F+ = W\{k + l)F+,i (13) 

G,k{enV) = -5!i±lzll!|p ; Fo+(£„p) = - (14) 

and we introduced renormalized frequency e and gap A as: 

?n = rjken] A = T^fcA; rik = l+ (15) 

lei + |A|2 



similar to those appearing for superconductors with impurities . 

From (p!2D-(p!5|) it is easy to obtain the system of recursion relations for real and imaginary 
parts of normal Green's function, as well as for anomalous Green's function: 



ImGk 
ReGk 



e — ImG 



{e - ImGf + {{-Ifip + ReGf + | A + F\^ 

{-l)% + ReG 

{e - ImGY + ((-l)^ep + ReGf + | A + F\^ 



F+ = = ^—^ = = ^ (16) 

^ (e-/mG)2 + ((-l)% + /?eG)2 + |A + F|2 ^ ^ 
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Let us introduce the following notations: 

ImGk = -BnJk; ReGk = -i-l)%Rk; = A* U (17) 

Then we can see that recursion relations for and fk just coincide, so that = f^. Finally 
we obtain the following system of recursion relations for and Rk- 

•Jk 



{el + A2)(ryfe + WH{k + 1) J,+i)2 + ^2(1 + W^v{k + l)Rk+iY 

o ^ 1 + + l)Rk+i . . 

' {el + /\'^){r,k + W^v{k + l)Jk+iY + il{l + W^v{k + l)Ru+iY ^ ' 

Then the normal and anomalous Green's functions of a superconductor are determined 
through Rq and Jq: 

ImG = -EnJo; ReG = -^pRo; = A* Jq (19) 

and represent the complete sum of perturbation series for electron in a superconductor 
interacting with AFM short - range order fluctuations. 



IV. CRITICAL TEMPERATURE AND TEMPERATURE DEPENDENCE OF 

THE GAP. 

Energy gap of a superconductor is determined by the equation: 

A{p) = -Tj2J2 Vscip, p')F{enp') (20) 

p' Sn 

On the fiat parts of the Fermi surface the anomalous Green's function is defined by (|19|) and 
recursion relations (piSl) . On the rest ("cold" part) of the Fermi surface the scattering on 
AFM fluctuations is absent (in our model), so that the anomalous Green's there function is 
given by (p^Ol). As a result, for the case of s-wave pairing, with the account of (H), the gap 
equation (|20|) takes the form: 



th 



l = X[c,Tj:lj(Me,.() + (l-&)l "(^^] (21) 

where A = VNq{0) is dimensionless coupling constant of pairing interaction (A^o(O) - free 
- electron density of states at the Fermi level ), a = Aa/n, where a is the angular size of 
a fiat patch on the Fermi surface (cf. Fig.l). In numerical calculations below we, rather 
arbitrarily, choose a = 2/3, i.e. a = tt/6, which is close e.g. to the data of Ref. |T2 . 



In case of d-wave pairing we have to take into account the angular dependence of the 
gap (^, so that Eq. (pO]) becomes: 




e^(0)E/ d^MenO+ d^e\^) f (22) 
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In Fig. 4 we show the calculated (from Eq. (0)) temperature dependences of the gap for the 
case of s-wave pairing and for different values of correlation length (parameter n = vfC,~^) 
of fluctuations. For the case of d-wave pairing results are qualitatively similar. 

Equation for superconducting critical temperature immediately follows from (pl|), (^2]) 
as A ^ 0. In this case Jo(^ 0) is independent of (p and is the same both for s and d-wave 
pairing. Then Tc - equation takes the following form: 

«e//Te^ / diUe^i; A ^ 0) + (1 - / di-^ (23) 

where an "effective" size of fiat patches on the Fermi surface is defined as: 

a = I " . ^ s-wave pairing) .^4) 

1 a + ^ sin(7ra) ( d-wave pairing) ' ^ ' 

Calculated dependences of Tc on the width of the pseudogap W and correlation length 
(parameter k = vf^~^) are shown in Fig.5 (T^o - transition temperature in the absence of 
pseudogap). 

The general qualitative conclusion is the same as in Refs. [|T0],^: pseudogap suppresses 



superconductivity due to a partial "dielectrization" of electronic spectrum on "hot" parts of 
the Fermi surface. This suppression effect is maximal for k = (infinite correlation length of 
AFM fluctuations) and weakens as correlation length becomes shorter. These results 

are in full accordance with experimental phase diagram of HTSC - cuprates. 

Let us stress once again that all our results are valid in case of self - averaging super- 



conducting order parameter (mean - field approach of Ref. [11|), which is valid for not very 
large correlation lengths ^ < where ^0 ^ is superconducting coherence length (the size 
of Cooper pairs at T = 0). For ^ ^ ^0 important effects due to non self - averaging gap 
fluctuations appear, leading e.g to characteristic "tails" in temperature dependence of the 



average gap for Tc < T < Tco |11 



V. COOPER INSTABILITY. RECURRENCE PROCEDURE FOR THE VERTEX 

PART. 

It is well known that the critical temperature can also be determined from Cooper in- 
stability of the normal phase: 

l-rx(0,0) = (25) 

where the generalized Cooper susceptibility is defined by the graph shown in Fig. 6. Here 
we have to calculate the "triangular" vertex part accounting for interaction with AFM 
fluctuations. For the similar one - dimensional problem (and for real frequencies T = 0) the 
appropriate recurrence procedure was formulated in Ref. p9[. For our two - dimensional 
model this procedure was used in Ref. to calculate optical conductivity. Generalization 



to Matsubara frequencies is rather straightforward. Below, for definiteness, we assume 
e„ > 0. Then we obtain: 
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{2ikn 1 
2ten - {-l)''VFq - W^v{k + l){Gk+i - Gfc+i) J 

T{en, -Sn, q) = ro(£„, -En, q) (26) 

where = Gk{£nP + q) and Gk = Gk{—£n, p) are calculated from (|]). 

To find Tc we have to know the vertex part at q = 0. Then Gk = Gl and the vertex Tk 
becomes real, which considerably simplifies (^). Using notations similar to (|T7|), from 
and (|26|) we get: 

while for i?^ and Jk we have recursion relations given by (|18|) with A = 0. 

There exists the following exact relation similar to the Ward identity (the proof of this 
relation will be given below): 

G(£„p)G(-£„p)r(£„, -Sn, 0) = (eJi?o(^nep) + elj'o{en^p))ro{en, -Sn, 0) = MSn^p) = J-!!!^^5^ 

(28) 

Numerical analysis fully confirms this relation, demonstrating the self - consistency of our 
recursion relations for one - particle Green's function and vertex part Q As Jo{A 0) 
coincides with Jq in normal phase, the relation (^) leads to - equation obtained from 
Cooper instability (^51): 

th^ 1 



1 = A ae//Te$^ r da^lRliSn^p) + £^ Jq' (^^.^p) )ro (^n, Sn, 0) + (1 - tte//) 

J —LUr JO 



e J 

(29) 

being the same as Eq. (^Bf), obtained by linearization of the gap equation, despite seemingly 
different recursion procedures used to obtain these equations. 

VI. GINZBURG LANDAU EXPANSION. 



In Ref. we derived the Ginzburg - Landau expansion in exactly solvable model of the 
pseudogap with infinite correlation length of AFM fluctuations. Here we generalize these 
results to the case of finite correlation lengths. 

Let us write the Ginzburg - Landau expansion for the difference of free energies of 
superconducting and normal state in the following form: 



^Note that an analytic proof of this relation from direct comparison of recursion procedures for 
Green's function and vertex part is non obvious, to say the least. 
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F,-Fn = AlAqp + g2c|Aq|2 + ||Aq|^ (30) 

where Aq - is the amphtude of the Fourier - component of the order parameter 

A(0,q) = A,e(0). (31) 



Now (|30D is defined by diagrams of loop - expansion for the free energy of an electron in 
the field of fluctuations of the order parameter with small wave - vector q [|T0[] . 
We express the coefficients of Ginzburg - Landau expansion as: 

A = AoKa; C = CoKc; B = B^Kb, (32) 

where Aq, Cq and Bq denote the standard expressions for these coefficients in the case of 
isotropic s-wave pairing: 

Ao = NM^: ^0 = iVo(0)|^|; = iVo(O)^, (33) 

Then all the anomalies of the model under consideration, connected with the appearance of 
the pseudogap, are contained in dimensionless coefficients Ka, Kq and Kb- In the absence 
of pseudogap all these coefficients are equal to 1, only in case of ci-wave pairing we have 
Kb = 3/2. Thus for d- wave pairing we shall appropriately normalize Kb, giving numerical 
results for Kb = 2/3Kb- 

Consider the generalized Cooper susceptibility shown in Fig. 6. 

x(qO;T) = -T ^ ^ G(£„,p + q)G(-£„p)e2(0)r(£„, q) (34) 

£n P 



Using (|2^) we can express the coefficients Ka and Kq as: 

X(q0;r)-x(00;r, 



K, 



A 







aeff-^h Yl r diUsnO-T, E r «(^„0l + l-«e// (35) 



ir, = lim.^^^°'^^)-^(°°'^^) 



2^ / d^JoiSnO 

e„=7rTc(2n+l) ''^'^ 



eff^ 



1 . 1 



-22 d^G{en,^ + -VFq)T{en,-en,q)G{-en,^ - -vpq)} + I - aeff (36) 

£=7rTe(2n+l) ^ ^ J 

Situation with coefficient B, in general case, is more complicated. Important simplifications 
appear if we limit ourselves in the order of |Aq|^, as is usually done, by considering only 
the case of g = 0. Then the coefficient B can be determined directly from the anomalous 
Green's function F, for which we already have the recurrence procedure (pISD , ([l9|) . Let us 
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consider the diagrammatic expansion of the anomalous Green's function, shown in Fig. 7(a). 
From this it becomes clear that: 

hm = G{enP)Gi-enP) + ... = G(£„p)G(-£„p)r(£„, -£„, 0) (37) 

which, by the way, immediately proves (|^) (use also ([19|) ). Thus, for the two - particle 
loop x(0, 0) we get: 

X(0, 0) = T E E li^^o = T E E Jo( A = 0) (38) 

P En P S„ 

For the "four - tail" diagram of Fig. 7(b), determining the coefficient B, in the same way we 
obtain: 

- TE E l-m ^ - '""r° ^ = E E l-n ^ °' (39) 

P en II P En II 



where Jo (A) is defined by the recursion procedure ([18|) . Finally, for the dimensionless 
coefficient Kb we have: 

Stt^TJ /-e Jo(A = 0)-Jo(A) 

where 

I a ( s-wave pairing) 

«B = i ~ , 4 • / ~A , 1 wo ) ^ • • ( • (41) 

a + 3^ sm(7raj + sm(27ra) ( d-wave pairmg) 

These expressions allow direct numerical calculations of the coefficients K^, Kc, Kb- In 
Fig. 8, for example, we present the calculated dependence of Kq on the width of the pseu- 
dogap W and correlation length of AFM fluctuations (parameter k, = vp/xi'"^). The appro- 
priate dependences of Ka and Kb are qualitatively quite similar. In particular, for k = 



we have just Kb = Kc [|10 



VII. PHYSICAL PROPERTIES OF SUPERCONDUCTORS WITH PSEUDOGAP. 

Ginzburg - Landau expansion deflnes two characteristic lengths of superconductors: the 
coherence length and penetration depth of magnetic fleld. 

The coherence length for given temperature ^(T) determines the characteristic scale of 
inhomogeneities of the order parameter A: 

eiT) = (42) 



In the absence of the pseudogap: 



escsiT) = (43) 
Cbcs{T) ^ 0.74 ^° , (44) 
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where — O.lSvp /Tc. In our model: 



e{T) 



e 



BCS 



(T) 



Ka 



(45) 



The dependences of ^'^{T)/^'j^Qg{T) on the width of the pseudogap W and correlation length 
of fluctuations (parameter k) for the case d-wave pairing are shown in Fig. 9. Note that the 
changes of coherence length are relatively small. 

For the penetration depth of a superconductor without the pseudogap we have: 



Ao 



%/2 Ji-T/r; 



^BCs{T) 



where Aq = is penetration depth at T = 0. In general case: 



(46) 



Then for our model: 



A(T) 



B 

■327re2 AC' 



Xbcs{T) \KaKc 



(47) 



(48) 



Graphical dependences of penetration depth for the case of d-wave pairing are shown in 
Fig.lO. 

Near Tc the upper critical magnetic field Hc2 is defined via Ginzburg - Landau coefficients 

as: 



H, 



A 



c2 



(49) 



where (po — cir/e is magnetic fiux quantum. Then the slope of the upper critical field close 
to Tc is defined as: 



dH, 



c2 



dT 



247r0o Ka 



dHc2 
dT 



(50) 



, normalized to the slope 



Graphic dependences of the slope of the upper critical field 

at Tco, on the effective width of the pseudogap W and parameter of correlation length n for 
the case of rf-wave pairing are shown in Fig. 11. It is seen that the slope for large enough 
correlation lengths rapidly drops with the width of the pseudogap. However, for short enough 
correlation lengths we can observe even some weak growth of this parameter for small values 
of the pseudogap width. For fixed pseudogap width the slope of i?c2 significantly grows as 
correlation length becomes smaller. 

Finally, let us consider the discontinuity of specific heat at the transition point: 



n 



B \T-Tr 



(51) 
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where Cs, Cn are specific heats of superconducting and normal states, Vt - sample volume. 
For T = Tco (in the absence of pseudogap, W = Gi): 

^■'-^"^ Niof-^. (52) 



^ ' ' 7C(3) 

Then the normalized discontinuity of specific heat in our model can be expressed as: 

{Cs — Cn)TcO ^cO Kb 



(53) 



Appropriate dependences on effective width of the pseudogap W and parameter of correla- 
tion length K for the case of d-wave pairing are shown in Fig. 12. It is seen that specific heat 
discontinuity rapidly drops with the growth of the pseudogap width and grows as correlation 
length of AFM fluctuations becomes smaller. 

For s-wave superconductor the dependences of physical properties are quite similar, the 
only change is in larger scale of W for which the appropriate changes appear, correspond- 
ing to larger stability of isotropic superconductors to partial "dielectrization" of electronic 
spectrum due to pseudogap formation on "hot patches" of the Fermi surface |]IO| , p!T| . 

Among the physical characteristics, analyzed above, relatively detailed experimental data 
are available for specific heat discontinuity In complete agreement with our conclusions, 
specific heat discontinuity for Bi — 2212 rapidly drops as the system moves to the underdoped 
region, where the width of the pseudogap grows. According to Ref. the width of the 
pseudogap (our parameter 2W) changes from the values of the order of 700K for hole 
concentration p = 0.05 to the values of the order of Tc ~ lOOK near optimal concentration 
p = 0.16 and drops to zero for p = 0.19. A relation between the drop of specific heat 
discontinuity and the growth of the pseudogap width is clearly observed. Unfortunately, 
we do not know detailed enough data on concentration dependence of correlation length of 
fluctuations and appropriate dependences of physical characteristics of superconductors in 
pseudogap state. Qualitatively it is obvious that correlation length grows as system moves 
to the underdoped region, so that the drop of specific heat discontinuity is natural also from 
this point of view. 



VIII. CONCLUSION. 



In this paper we continued the study of anomalies of superconducting state in the frame- 
work of rather crude model of the pseudogap state of two - dimensional electronic system 
10| , p]T|| , which is, however, in qualitative agreement with a number of observed anomalies of 
electronic structure of underdoped HTSC - cuprates. In Refs. |]TO|,|n[ we considered rather 
unrealistic limit of infinite correlation length of AFM short - range order fluctuations, which 
allowed us to obtain an exact analytic solution. Here we generalized our model to realistic 
case of finite correlation lengths, taking into account, as in Refs. [0,|ll[], all diagrams of 
perturbation theory on electron interaction with fluctuations of short - range order. Our 
analysis was performed in a standard (mean - field in terms of Ref. [|1T|) approach, assum- 
ing self - averaging property of superconducting order - parameter over fluctuations of the 



random field of AFM fluctuations. In Ref. |TT| we have shown that this assumption is not 
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justified in the limit of ^ ^ oo. At the same time this assumption is apparently well justified 
for the case of ^ -C (where is the coherence length of a superconductor at T = 0, i.e. 
the size of Cooper pairs). Thus, we are left with rather complicated task of the accounting 
of non self - averaging effects for ^ > C,o- We have already mentioned that in real HTSC 
systems in most cases ,^ ~ .^o, so that effects of non self - averaging superconducting gap, 
similar to those considered in Ref. may be very important, leading e.g. to characteristic 



"tails" in the temperature dependence of the average gap for T > Tc and the physical picture 



of superconducting "drops" of Ref. [11 



Another serious simplification of our model is the assumption of static (and Gaussian) 
nature of short - range order of fluctuations. This assumption may be justified only for high 
enough temperatures T ^ Usf (where Ugf - is characteristic scale of spin fluctuations) [^,0 . 
Thus, the use of static approximation in superconducting state for T < Tc is rather doubtful. 
However, we think that our simplified treatment allows us to describe most important effects 
of the changes of the electronic spectrum (due to pseudogap formation on "hot patches" of 
the Fermi surface) upon superconductivity. The account of spin dynamics inevitably requires 
to drop the simple phenomenology of BCS model and consider the microscopic nature of 
pairing interaction. It is doubtful that such a program can be realized in near future. In 
particular, the problem of summation of all perturbation theory diagrams for the interaction 
with dynamical spin fluctuations seems absolutely hopeless. 

This work was partially supported by grants 99-02-16285 of the Russian Foundation of 
Basic Research and REC-005 of CRDF, as well as by the projects 108-11(00-11) of the State 
Program "Statistical Physics" and 96-051 of the State Program on HTSC. 
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APPENDIX A: COORDINATE REPRESENTATION. NORMAL AND 
ANOMALOUS GREEN'S FUNCTIONS. 



Let us consider some technical details of derivation of recursion relations for Gorkov 
equations (|12]) - (pTSl) . It is sufficient to limit the analysis to consideration of two flat 
parts of the Fermi surface, orthogonal to Px - axis, which are connected by the scattering 
vector Q = {±2pp, 0). Then the problem becomes purely one - dimensional as the velocity 
projection Vy = and the linearized electronic spectrum ^p^^fpp = ±.vfPx does not depend 
on ^/-component of electronic momentum. For brevity in the following we just put Vp = 1. 

Calculations simplify in coordinate representation |2T| , considering the electron propaga- 
tion if the field of Gaussian AFM fiuctuations W{x) ^ W*{x) (incommensurate case) with 
correlator: 

< W*ix)W{x') >= iy2g-.k-x'| (A^) 

Then electron propagators corresponding to normal and anomalous Green's functions of a 
superconductor (p^, take the form: 



Mx) = r ^e^'^'^F.^Px) = , ^* e"V^^I-l (A2) 

/el + lAP 



where = 1 for right moving particles, and (73 = —1 for left moving particles. Scattering 
by fiuctuations transforms "right" particles into "left" and vice versa. From ([A2|) it is seen 



that the particle moving along the path of the length I produces the factor e~ v^^^'. 

During calculation of specific diagrams it is convenient ||2l| to change integration variables 
from coordinates of interaction vertices Xk to lengths of the paths Ik traversed by the particle 
between separate scatterings, fixing the total displacement x — x' . With interaction line 
connecting the vertices m and n on electronic line we have to associate the factor: 

W''exp{-K\xm - Xn\) = W\xp{-K\ J] (A3) 

k=ni 

Integration over all Ik is performed from to 00. 

It is seen that the finite correlation length of fiuctuations in a given diagram leads to some 
"damping" of the appropriate transition amplitude with distance traversed by an electron. 
The exact treatment of this effect is difficult, but in Ref. we used an obvious inequality: 

n— 1 n— 1 

exp(-K\ (-1)''4|) > exp[-K Ik) (A4) 

k=m k=m 



and replaced the exponential of (|A3|) by the exponential from the r.h.s. of (|A4|) . This 



is equivalent to the replacement of correlator of random fields ( [Al| ) by similar expression, 
where in the exponent we just replace the distance |x — x'l by the total length of the path 
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traversed by the particle between scattering acts at x and x' . According to ( |A4|) this proce- 
dure somehow overestimates the effect of damping k in each diagram of perturbation theory. 
After such a replacement the diagrams of all orders are easily calculated and precisely repro- 
duce the Ansatz of (|^) for the normal phase ||2l|. We have already mentioned that results 



obtained in this way, e.g. for the density of states, are in good quantitative agreement with 
exact numerical simulation of the one - dimensional problem [^,^ , which gives additional 
support for our approximation strengthening qualitative estimates of Ref . |^ . 

Let us use the same approximation during the analysis of diagrams of perturbation theory 
in superconducting phase, which are built upon propagators (|A2|) . In this case the role of 



interaction reduces just to the appearance of additional factor of e in each normal or 



anomalous Green's function ( |A2D , surrounded by the given interaction line, or (which is the 



same) to the addition of k to y + | in the exponential of each Green's function. Making 
transformation back to the momentum representation it is easily seen that the contribution 
of an arbitrary diagram of the higher order of perturbation theory is determined by the 
product of the appropriate number of normal and anomalous Green's functions of the form: 

GM = ^^"^'f ' ; F^,{p)= ff ^^ ■ (A5) 



where 5^ = ye^ + |Ap + /c/t, while A; - is the number of interaction lines, surrounding the 
given Green's function. The factor of (—1)^ is due to scattering transforming "right" particles 
into "left" and vice versa. Introducing the renormalized frequency and gap as in (|15|) , we 
can see that ( |A5| ) reduces to the standard form which completes the justification of 
our recursion procedure (|12]), (pISl). 



16 



Figure Captions: 



Fig.l. Fermi surface of two - dimensional system. "Hot patclies" are sliown by tliick 
lines of the width ~ 

Fig. 2. Diagrammatic representation of recursion relation for one - particle Green's func- 
tion. 

Fig. 3. Diagrammatic representation of recursion relations for Gorkov equations. 

Fig. 4. Temperature dependence of superconducting energy gap for the case of s-wave 
pairing and different values of correlation length (parameter k = vf^~^) of AFM fluctuations, 
calculated for A = 0.4, ^ = 3: 

f=0(l); 1.0(2); 10.0(3). 

Dashed line — A(T) in the absence of the pseudogap. 

Fig. 5. Dependence of superconducting transition temperature on the width of the pseu- 
dogap W and correlation length of AFM fluctuations (parameter k = vf^~^)- 
f =0.1(1); 1.0(2); 10.0(3). 
Dashed line — k = |T0|. 

At the insert: dependence of Tc on k for ^ = 5. 

Fig. 6. Diagram for the generalized Cooper susceptibility. 

Fig. 7. (a) - Diagram series for anomalous Green's function, dashed lines - AFM fluctu- 
ations, (b) - Diagram determining Kb- 

Fig. 8. Dependence of Kc on the width of the pseudogap W and correlation length of 
AFM fluctuations (parameter k = vfC,~^)'- 

f =0.1(1); 1.0(2); 10.0(3). 

Dashed line — n = |10[ . 



At the insert: dependence of Kc on k for 7^ = 5. 

-t cO 

Fig. 9. Dependence of supercoducting coherence length on the width of the pseudogap 
W and correlation length of AFM fluctuations (parameter n = vfC,~^)'- 
§=0.l{iy, 1.0(2); 10.0(3). 
Dashed line — k = [|T0[ . 



At the insert: dependence of coherence length on k for ^ = 5. 
Fig. 10. Dependence of penetration depth on the width of the pseudogap W and correla- 
tion length of AFM fluctuations (parameter k = vfC,~^)'- 
f =0.1(1); 1.0(2); 10.0(3). 

At the insert: dependence of penetration length on k, for ^ = 5. 

J- cO 

Fig. 11. Dependence of the slope of the upper critical field on the width of the pseudogap 
W and correlation length of AFM fluctuations (parameter k = vf^~^)'- 
#=0.1(1); 1.0(2); 10.0(3). 
Dashed line — k = [|10[ . 



At the insert: dependence of the slope of Hc2 on k for ^ = 5. 
Fig. 12. Dependence of specific heat discontinuity on the width of the pseudogap W and 
correlation length of AFM fluctuations (parameter n = vfC,~^)'- 
f =0.1(1); 1.0(2); 10.0(3). 
Dashed line — k = |jlO[ . 



At the insert: dependence of specific heat discontinuity on k for ^ = 5. 
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